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Abstract 

Let q > 1 be an integer. Given M samples of a smooth function of q variables, 27r-periodic 
in each variable, we consider the problem of constructing a q-variate trigonometric polynomial of 
spherical degree 0(M 1 ^ q ) which interpolates the given data, remains bounded (independent of M) on 
[— 7r,7r] 9 , and converges to the function at an optimal rate on the set where the data becomes dense. 
We prove that the solution of an appropriate optimization problem leads to such an interpolant. 
Numerical examples are given to demonstrate that this procedure overcomes the Runge phenomenon 
when interpolation at equidistant nodes on [—1, 1] is constructed, and also provides a respectable 
approximation for bivariate grid data, which does not become dense on the whole domain. 



1 Introduction 

Interpolation at equidistant nodes on the unit interval [—1, 1] is a very classical problem. In the first course 
in numerical analysis, one learns of the Newton divided difference algorithm to find such an interpolant, 
and the corresponding error formula. The Runge example, x (x 2 + 25) _1 , shows that the sequence 
of these interpolants need not converge even if the target function is analytic on [—1,1]. In general, 
Faber's theorem [13l Theorem 2, p. 27] states that for any interpolation matrix on [—1, 1], there exists a 
continuous function on [— 1 , 1] such that the corresponding polynomials of interpolation to this function 
do not converge. 

The situation changes drastically if one allows the degree of the interpolatory polynomial to be greater 
than the minimal required. Thus, the following Theorem 1 1.1 1 is a simple consequence of [17| Theorem 2.7, 
p. 52]. For the purpose of this exposition, we denote the class of all algebraic polynomials of degree at 
most m by LI m , and define ||/||oo.[-i.i] : = su Pte[-i i] I /Ml- We note that for n equidistant nodes on 
[— 1, 1], the quantity d n in the following theorem satisfies d n > 2/n. 

Theorem 1.1 Let Xk, n = cosOk,n G [ — 1, 1] be an arbitrary system of nodes (0 < d\_ n < • • • < 9 n _ n < it) 
and let 

d n := mm {9 k +i, n - 0k, n )- 

l<fc<n-l 

Then for any e > 0, there exist linear polynomial operators P n on C[— 1, 1] with the following properties: 
(a) If m = [7r(l + e)/d n \ then P n (P) = P for all P G U mi (b) for f G C[-l, 1], P n (f) £ Il N where 
N = (7r/d„ + l)(l + 3e), (c) P(f, x k , n ) = f(x k>n ) for k = 1, • • • , n, and (d) 

11/ - P«(/)|U,[-i,i] < c inf ||/ - P|U,[-i,i]. (1.1) 
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In many engineering applications, one has to find a good approximation to an unknown multivariate 
target function which also interpolates the function at certain points, sometimes called landmarks. For 
example, in the problem of image registration, we are given a set of locations Xj G [— 1, l] 2 in the first 
image and a corresponding set of points yj G [— 1, l] 2 in the second image. The idea is that the location 
xj in the first image is the "same" as the location y 3 in the second image. We then hope to find a map 
g : [—1, l] 2 — » M 2 such that g{xj) — yj, and such that g satisfies some smoothness conditions. There are 
at least two reasons for insisting on interpolatory approximation in this situation. First, the locations 
might have been chosen at great costs, including human efforts. Second, if the registration is being done 
many times over a sequence of images (for example when we stitch together video frames to form a large 
image), then a non- interpolatory approximation will cause a drift between the first image and the last 
image in the sequence. 

It is interesting to note that polynomial interpolation in multivariate setting has a totally different 
flavor than in the univariate setting; for example, even if one has exactly as many points as the dimension 
of the polynomial space involved, there might not exist an interpolant from that space. Even if an 
interpolant exists, the error bounds for approximation depend heavily on the geometry of the points. In 
[8] , we proved that an analogue of Theorem If .11 holds in practically any setting where the so called direct 
theorem of approximation holds, provided we drop the requirement of linearity. In particular, we proved 
analogous results in the multivariate setting. However, the results in [5] are not constructive, and do not 
yield linear operators. 

The purpose of this paper is to develop algorithms to achieve near best polynomial approximations 
to smooth multivariate functions, which satisfy interpolatory constraints. Our constructions will work 
without requiring any specific locations for the points where the target function is evaluated. We refer to 
such data as scattered data. We do not require that the data become dense on the whole cube. In turn, 
our approximations may not converge on the whole cube. However, they will converge at the limit points 
of the data, and we will estimate the rate of convergence. 

To motivate our construction, we revert to the univariate case of Theorem ll.il We recall that there is 
a one to one correspondence between functions on [—1, 1] and even, 27r-periodic function on R, given by 
= /(cos0). Moreover, || || oo, [— := sup ee[ _ W 7r] |/°(0)| = ||/||oo,[-i,i]- Let r > 1 be an integer, 
and f° be r times continously differentiable on [— 7r,7r]. In this discussion, we will write P n in place of 
P n (/). In view of a theorem of Czipser and Freud [2], the estimate implies that ||-P„ ^ ||oo,[— tt.tt] < 
cll/ ^ ||oo. [-7t.it]- Therefore, the minimization problem "minimize ||P°( r ) ||oo,[-7r,7r] over all P G Iljv, 
subject to the constraints P(xk. n ) = f{xk,n), k = 1, •••,n" has a solution P* with the right bounds 
on the r-th derivative of its periodic version. In view of the Arzela-Ascoli theorem, this implies that 
any subsequence of the sequence {P^} has a uniformly convergent subsequence. If xo is a limit point 
of a subsequence {xk,n}neA, then it is not difficult to deduce using the interpolatory conditions that 
lim n6 A,„_ s . 00 P* (xq) = f{xo). In this paper, we will extend these ideas to the multivariate periodic setting. 
Instead of describing the smoothness of the functions in terms of derivatives, we will consider Sobolev 
classes. We will also consider minimization in arbitrary L p norms; the L 1 norm being of recent interest 
from the point of view of compressed sensing. Some technical details, involving a construction of quasi- 
interpolatory polynomial operators, are required to prove the rate of convergence of our constructions. 
However, the bulk of the technical details is in the proof of the feasibility of the optimization problem. 
We will use Theorem 2.1 in [5] with the appropriate Sobolev spaces, and will need to prove the analogue 
of Theorem 3.2 in [8] also with approximation in Sobolev spaces rather than the space of continuous 
functions as in that theorem. Our main tool is the construction of a multivariate analogue of trigonometric 
polynomial frames constructed in |11[ 112) . 

We state our main results Section [3J and illustrate them numerically in Section [3] The proofs of the 
results are given in Section [SJ following some preparation of a technical nature in Section At the first 
reading, it might help to skip this section, referring back to the various statements there on an as needed 
basis. 

We would like to thank Karthik Raghuram for carrying out all the numerical experiments. 
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2 Main results 



In the sequel, q > 1 will denote a fixed integer, and we will think of 2-7r-periodic functions on R 9 
as functions on [— 7r, 7r] 9 , tacitly identified with the q dimensional torus. Analogous to the univariate 
case, any function / : [— 1, l] q — > R, corresponds uniquely to the 27r-periodic function f° on R 9 by the 
correspondence 

/ (#!,••• A) =/(cOS0i.---,COS0,). 

The symbol || o || will denote the Euclidean norm of a vector in R 9 . Let H 9 denote the class of all 
trigonometric polynomials in q variables with spherical order at most n; i.e., 



Yl akexp(ik-(o)) : a k G C}. 



keZ 3 , ||k||<n 

Here, we find it convenient to use the same notation even if n is not an integer. It is not difficult to 
see that multivariate algebraic polynomials on [— 1, l] 9 correspond to the trigonometric polynomials of 
the same order which are symmetric in each of the variables. Therefore, in this paper, we are interested 
mainly in the interpolation of multivariate periodic functions; the results can also be applied trivially to 
the interpolation of functions on [— 1, l] 9 , with suitable smoothness conditions defined in terms of the 
corresponding periodic function. 

If 1 < p < oo, K C [— 7T, 7r] 9 and / : K — >• C are Lebesgue measurable, we write 

{( / |/(x)| p dxl , if 1< p< oo, 
Uc IA n J ~ P (2.1) 

ess sup |/(x)|, ifp = oo. 

The symbol L P (K) denotes the class of all Lebesgue measurable functions / for which ||/|| P! _k- < oo, with 
the usual convention that two functions are considered equal if they are equal almost everywhere. If 
K = [— 7T, 7r] 9 , we will omit its mention from the notations. If 1 < p < oo, we will write p' :— p/(j> — 1), 
and extend this notation to p = l,oo by setting 1' = oo, oo' = 1. If / € L 1 , the Fourier coefficients of / 
are defined by 

/(k) := -1- f /(x) expHk ■ x)dx, k G Z 9 . (2.2) 
If / e L p , then its degree of approximation from H 9 is defined by 

EnM) := #4 ll/_T|lp - 

If s e R, 1 < p < oo, the Sobolev class Wf consists of all / e LP for which there exists /' s ^ g L p such 
that 

/M(k) = (||k|| 2 + ir/ 2 /(k), kez 9 . 



We define 



\f\\wS :=II/ (S) II P , (2-3) 



and note that Wf is a Banach space. We observe that if A is the Laplacian operator on R 9 , and s is 
an even, positive integer, then /^ s ' = (A + I) s / 2 f, where / is the indentity operator. In particular, in 
this case, the operator / H> is a surface derivative operator on the torus identified with [— tt, n] q . 
An important property of the spaces Wf is given in the following proposition, which will be proved 
in Section l4~2l Here, and in the rest of this paper, the symbols c, ci,--- will denote generic positive 
constants, depending on such fixed parameters of the problem as p, s, q, etc. and other quantities 
explicitly indicated, but their value may different at different occurrences, even within a single formula. 
The notation A ~ B means that c%A < B < c^A. 

The following proposition, to be proved in Section 14.21 gives an integral representation of functions 
in Wf. 
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Proposition 2.1 Let 1 < p < oo, s > q/p. Then there exists a function K s £ L p such that 

K s (k) = (||kf + l)- s / 2 , keff. (2.4) 
// / £ Wf, then for almost all x <E [— tt, ir) q , 

f( x ) = J—[ Ksix - y )fM {y)dy= J_r K s (y-x)f^(y)dy. (2.5) 

In particular, f is almost everywhere equal to a continuous function. Denoting this continuous function 
again by f , we have for any < s' < s — q/p, 

£ 2 n,oo(/) < c2-"(-«/rt||/|| w? , ll/IU < c||/|| w - < c\\f\\ wi . (2.6) 

We remark that in the case p = 2, one can take the following approach for interpolation of functions 
in Wg, s > q/2. The Golomb-Weinberger variation principle [3] can be used to show that the solution of 
the minimization problem 

minimize : g{y ,n) = f(yj,n), j = h" -, M n} (2.7) 

has a solution in the span of -K~2s(° ~ Yj.n), an d therefore, can be found by solving an appropriate system 
of linear equations. The stability of this system as well as the error bounds can be estimated using known 
techniques from the theory of radial basis functions, for example, [14] (See Theorem 14.51 below) . However, 
we are interested in finding polynomial interpolants for functions in Wf for s > q/p, without requiring 
p = 2. 

As customary in the theory of interpolation, let Y be the interpolation matrix whose n-th row Y n 
contains M n vectors {yj, n }f=i- Our theorems will depend upon two quantities, defined in (|2.8[) below, 
that measure the density of these points as well as their rareness. First, if C C [—n, ir] q and x £ [—tt, ir] q , 
we define 

dist(C,x) := inf ||x — y||. 

Further, if K C [— tt, n] q , we define the mesh norm d(C,K) (respectively, separation radius 77c) of C by 
S(C,K) := supdist(C,x), r] C ■= (1/2) inf ||x-y||. (2.8) 

We will simplify our notation, and write 8 n {K) for 5(Y n , K) and r\ n := n(Y n ). 

Our first theorem, to be proved in Section [5l shows the feasibility of a procedure for finding interpo- 
lator trigonometric polynomials in H*. 

Theorem 2.1 Let 1 < p < 00, s > q/p, Y be as above. 

(a) There exists an integer N* with N* ~ i]^ 1 and a mapping P : Wf — > H^* such that for every f £ Wf , 

P(/, Yj.n) = f(y j<n ), j = 1, • • • , M n , (2.9) 

and 

\\f-P(f)\\ wF <cM{\\f-T\\ m : T£U%,}. (2.10) 

(b) We consider the minimization problem 

!i/ P 
(2.11) 



minimize 



where the minimum is over all P £ M q Nt , such that P{yj t n) = f{yj.n), j = 1,' • • ■> M n , and an ap- 
propriate interpretation is understood in the case p = 00. There exists a solution of this problem, 
P* = ¥* n (p,Y n ,f) £ W q N ,, such that \\P* n \\ W p < c\\f\\ wF . 
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In practice, it seems that we can take N* = 4r/~ 1 if p = 2 and q = 1. We note that the problem (|2.11|) 
has a unique solution if p = 2, and the corresponding operator P* is linear in /. 

The next theorem, to be proved in Section [SJ examines the convergence properties of the sequence 

{K}- 

Theorem 2.2 Let 1 < p < oo, s > q/p, f G Wf , N* and P* be found as in Theorem 0OI 

(a) If A is a subsequence of positive integers, xo G [n,Tt] q , and 

Jim dist(r„,x ) = 0, (2.12) 

ngA 

then 

Jim P;(xo) =/(x ). 

ngA 

(b) There exists a constantly = "f(p, q, s) (independent of n) with the following property. i/xo G [— n,Tt] q , 
and <5 n ([xo — <5, xo + S] q ) < 7<5, i/ien 

11/ - P;iU,[x - 5)X0+5 ], < c5 s -^\\f\\ wi . (2.13) 

The proof of Theorem 12 . 1 1 occupies a major part of this paper. We will use an abstract result from [5], 
quoted here as Lemma l5.ll To use this result, we need first to approximate carefully an arbitrary element 
of the span of {K s (o — y) : y G Yn*} for a suitable value of N* by trigonometric polynomials in Hat*; 
indeed, N* will be determined so that this approximation works. In turn, this involves an estimation of 
the coefficients of this element in terms of the norms of this element, as well as a good approximation 
bound on K 8 . In preparation, in Section 14.11 we introduce certain localized kenels and operators, and 
prove a number of technical results concerning these. These enable us after some further preparation to 
prove Proposition 12.11 and study some further properties of the kernel K s in Section 14.21 The proof of 
Theorem 12.2( a). as expected, is a compactness argument. We also need to estimate the discrete norm 
used in (|2.11|) by the corresponding continuous norm. The necessary facts are stated in Lemmas 14.41 and 
14.31 The proof of Theorem |2.2f b) is quite simple in the case when q = 1, p = oo, and s in an integer: If 
there are s elements of Y n in K = [xq — 5, xq + 8] , we take the Lagrange interpolatory polynomial L for 
/ (and hence, Pat*) at these points. The elementary Newton error formula for interpolation yields 

11/ - PiV* IU.JC < 11/ - L\\oo,K + ||Ptf. ~ LIU.JC < c5 S ||/ (s) ||oo- 

The Newton formula does not hold in the multivariate case, and no similarly clean estimates are possible 
independently of the geometry of the points in question. Therefore, we use a result from [9] , quoted as 
Proposition 14.51 below, to construct an analogue of L, which is not interpolatory, but utilizes only the 
values f(yj, n ) = Pjv*(yj,n)- To take care of the technicalities of nonintcgcr s and LP norms other than 
p = oo, we use the direct and converse theorems of approximation theory. Although these results are 
folklore, we could not find them in the literature in the form which we needed. Therefore, in Section [4.31 
we review the results in the form in which we found them, and reconcile them to our needs. We also 
describe the construction of the algebraic polynomial approximation. 



3 Numerical experiments 

In this section, we will present numerical experiments that demonstrate the behavior of the method over 
a wide variety of situations, some of which do not satisfy the assumptions made in this paper. In the case 
when p = 2, the optimization problem (|2.11j) has a numerically effective closed-form solution. In this case, 
the problem is formulated easier directly in terms of the coefficients of the trigonometric polynomials: 

Find argmin ^ |a k | 2 (l + ||k|| 2 ) s , (3.1) 
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subject to the constraints 

^2 a k exp(ik • y J: „) = f(yj, n ), j = l,---, M n . 

||k||<iV* 

Therefore, we assume in this section that p — 2, and refer the interpolant resulting as a solution of 
this problem as a minimum Sobolev norm (MNS) interpolant. In our computations below, we actually 
consider rectangular sums rather than the spherical sums as in (|3.ip . The term MNS interpolant will be 
used for all such minor variations. 

We first consider the classical Runge phenomenon by interpolating the function f(x) = (1 + lOOcc 2 ) -1 
at equi-spaced points {—1 + 2j/(n + l)}"=i on the interval [—1, 1]. In Table Q] we show the results of 
our numerical experiments. The first column of the table shows the number of data points n used for 
interpolation. To avoid any special structure among the points, we chose the values of n as indicated, so 
as to be essentially (but not exactly) doubling from step to step. In all cases, we chose the order of the 
intcrpolatory polynomial to be 2n, and computed the maximum error by sampling the MSN interpolant 
at 3n equi-spaced points. Columns 2-8 show the maximum interpolation error with different values of 
s. The maximum error decreases more rapidly with increasing s, but there are diminishing returns for 
higher values of s due to increasing condition numbers and the concomitant loss of numerical accuracy. 
To minimize this loss, we used a special algorithm that combines an LU factorization along with the 
traditional LQ factorization for solving the minimum norm problem [T]. Our computations show clearly 
that the interpolants converge; i.e., the Runge phenomenon has disappeared. For comparison we also 
show in the last column the approximation error from using a cubic-spline interpolant. 



n 


s = 1.5 


s = 2.5 


s = 3.5 


s = 4.5 


s = 5.5 


s = 6.5 


Spline 


31 


3.6212e-03 


3.1000e-03 


4.2114e-03 


5.1510c-03 


1.0535e-01 


1.0127e+00 


3.5710e-03 


61 


4.5758e-04 


1.0681c-04 


4.7994e-05 


3.0439e-05 


2.4317e-05 


2.0644c-04 


6.5167e-04 


121 


1.78446-04 


1.2509e-06 


8.4610e-08 


3.2494e-08 


5.1699e-08 


2.7949e-07 


4.1035e-05 


241 


6.7863e-05 


2.4084e-07 


6.8299e-09 


1.1385e-10 


2.5882e-09 


2.7632e-07 


2.3897e-06 


481 


2.5210e-05 


4.5175e-08 


6.4371e-10 


5.8196e-ll 


6.5261e-09 


7.2588e-06 


1.4739e-07 


961 


9.2203e-06 


8.3199e-09 


6.3179e-ll 


7.5275e-ll 


6.2423e-08 


7.3983e-04 


9.2473e-09 



Table 1: Columns 2-7: maximum error for MSN interpolant of Runge's function (1 + 100x 2 ) -1 on [—1, 1] 
with n equi-spaced points and different choices of s. The order of the interpolant was 2n in all cases. 
Column 8: maximum error with cubic spline interpolant. 

Next, we consider a two dimensional interpolation problem on a region inside the square [—1,1] x 
[— 1, 1]. The target function is given in polar coordinates by 

f(r,9) := |r-l/4| 1/8 |l-r| 4/5 sin(r(2cos6> + sin<9)). 

The function is singular on the circles of radii 0, 1/4 and 1. Furthermore the function does not satisfy the 
smoothness conditions of this paper. For the data points, we take those vertices of a square grid which lie 
in the indicated regions. If h is the length of each side of the squares in this grid, the target polynomial 
is a bivariate polynomial of coordinatewise degree [2/h\. In the following tables, n denotes the number 
of grid points which lie in the region in question, and m is the dimension of the space of interpolatory 
polynomials. 

In Table [2] we compute the maximum error of the interpolant in the annulus 1/2 < r < 3/4 using 
approximately An grid points. The results are shown in Table [5J This particular annulus is well removed 
from the singularities of /. Therefore it is pleasing to see that the MSN interpolant approximates the 
underlying function very accurately. We report the maximum error in the annuli 3/4 < r < 19/20 
and 1/4 < r < 3/10 for the same MSN interpolants in Tables [3] and 0] respectively. These annuli are 
significantly closer to the circles of radii r = 1 and r = 1/4 where the function is singular. Not surprisingly, 
the error is much larger here, but still usefully small. 

Next, for the same function /, we restricted the samples to the region {{r < 1/4} U {1 < r}} n 
{[—1,1] x [—1,1]}. Note that this region is essentially made up of 5 pieces. We used n equi-spaced 
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n 


m 


s = 1 


s = 2 


s = 3 


s = 4 


s = 5 


100 


484 


2.1979e-03 


1.4012e-03 


2.0732e-03 


2.1384e-03 


3.2363e-03 


352 


1764 


9.0618e-04 


1.0095e-03 


1.5181e-03 


2.7093e-03 


5.0088e-03 


1280 


6724 


2.4002e-04 


1.4093e-04 


2.0451e-04 


2.1890e-04 


3.4952e-04 


4924 


26244 


2.3078e-04 


1.0745e-05 


2.3017e-05 


3.4346e-05 


1.4726e-04 



Table 2: Maximum error of MSN interpolant in region 1/2 < r < 3/4. 



n 


m 


s = 1 


s = 2 


s = 3 


s = 4 


s = 5 


100 


484 


2.8778e-02 


4.6506e-03 


1.2650e-02 


3.3905e-02 


8.8795e-02 


352 


1764 


3.5019e-03 


6.8887e-04 


5.4859e-03 


2.5252e-02 


7.2012e-02 


1280 


6724 


1.4967e-03 


3.3349e-04 


7.7618e-04 


6.3416e-03 


2.9345e-02 


4924 


26244 


1.4569e-04 


1.8086e-04 


4.1918e-04 


9.2673e-04 


3.2014e-02 



Table 3: Maximum error of MSN interpolant in region 3/4 < r < 19/20. 



samples in the region and computed the MSN interpolant with m coefficients for different values of s. 
The maximum error in the region r < 1/5 is reported in Table and in the region 1.1 < r in Table [6] 

These experiments show that the proposed scheme can perform well even on difficult problems, es- 
pecially in two dimensions where traditional interpolation schemes require much more work to achieve 
comparable accuracy. The proposed method requires special algorithms to execute efficiently, which will 
be discussed elsewhere. The ideas presented here can also be generalized to handle noisy and redundant 
observations. These matters will also be reported elsewhere pQ. 



4 Technical preparation 

In this section, we present may technical results which are preparatory to the proof of the main results of 
Section [5] Our proof of Theorem 12.11 will require Theorem 14.21 and Theorem 14.41 Subsections 14. II and 14. 21 
are devoted to the proof of these. In Subscction l4.11 we introduce a localized kernel and the corresponding 
operator which will be used throughout this paper, and prove a number of results regarding these. In 
particular, we use these results in Subsection 14.21 to prove Proposition 12.11 and establish a few other 
facts related to the kernel K s . In Subsection I4.3[ we review some well known properties of multivariate 
trigonometric and algebraic polynomial approximation, which will be used in the proof of Theorem 12.21 



4.1 Localized kernels 

Let q > 1 be an integer. For t > 0, and h : [0, oo) — > R, we define formally 

%(h,x) := ^2 /i(||k||/t)exp(ik-x), xel«. (4.1) 

We set *oO,x) := 1 and %(h,x) := if t < 0. 

The following theorem summarizes the important localization estimate for the kernel where we 
use the notation 

Vf{u) = f(u)/u. 

Theorem 4.1 Let Q > (q + l)/2 be an integer, h : [0, oo) — > [0, oo) be a Q — 1 times continuously 
differentiable function supported on [0,1], with an absolutely continuous derivative h^~ x \ In addition, 
we assume that for some constants < a < b < oo, h(t) =0 if t > b, and h'(t) = if < t < a. With 
R = (q — l)/2 + Q, we have 

\M^x)\<4V Q Hi,[0,oo) . n *L . hflV xG[-ir,7r]«, t>0. (4.2) 

mm(l, (tHxIl)-") 
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n 


m 


s = 1 


s = 2 


s = 3 


s = 4 


s = 5 


100 


484 


3.1104e-02 


2.9877e-02 


2.9187e-02 


2.9011e-02 


2.8334e-02 


352 


1764 


6.1838e-02 


5.2437e-02 


4.8763e-02 


4.6610e-02 


4.5226e-02 


1280 


6724 


1.2503e-02 


8.0592e-03 


7.8406e-03 


8.6545e-03 


9.2571e-03 


4924 


26244 


2.0597e-02 


1.4048e-02 


1.0375e-02 


8.6748e-03 


8.6535e-03 



Table 4: Maximum error of MSN interpolant in region 1/4 < r < 3/10. 



n 


m 


s = 1 


s = 2 


s = 3 


s = 4 


s = 5 


21 


484 


7.7266e-03 


1.4501e-02 


1.7311e-02 


1.9202c-02 


2.0064e-02 


89 


1764 


2.1075e-03 


1.8032e-03 


2.0722e-03 


2.1088e-03 


2.1234e-03 


401 


6724 


2.2433e-03 


1.2778e-03 


7.6602e-04 


5.5885e-04 


4.8267e-04 


1637 


26244 


8.9254e-04 


9.2024e-04 


5.7556e-04 


2.9398e-04 


1.5738e-04 



Further, 

and for I < p < oo . 



Table 5: Maximum error of MSN interpolant in region r < 1/5. 



max \%(h, X )\=^ t (h,0)<ct"\\V^h\\ lA0)Oo) , t > 0, (4.3) 



||^(/i,o)|| p <^/f ||X>O/i|| 1i[0jOo) i>0. (4.4) 
Here, the constants denoted by c may depend upon a, b, q, and Q only. 

In order to prove this theorem, we recall that the Bessel function J a can be defined for a > —1/2, 
t > by ([HI Formula (1.71.6)]) 

/ (t) - ^/ 2 )" f 1 e ««d _ u 2 ) a - 1/2 du 

M > ~ r((2a + i)/2)r(i/2)l 1 e (1 uj 



(t/2)< 



r((2a + l)/2)r(l/2) 

It is customary to define 



l 
l 

e~ 



ltu (l~u 2 ) a ~ 1/2 du. (4.5) 



1/2) 

For / 6 L 1 (R ? ), we define its inverse Fourier transform by 



J -i/2 (*) = f,(iV cos *> f > °- ( 4 - 6 ) 



/(x) = (2tt)-« / /(y) cxp(zy • x)dy, x € f. (4.7) 



Lemma 4.1 (a) Let h (x) = h(\\x.\\), x e M«. Tften 

x||( 2 '^)/ 2 /■ 
(2tt)«/ 2 A 



~ llxll ( 2 -«)/ 2 r°° 

feo(x) = " " . /2 / M*)./( 9 - 2 )/2(s||x||) S «/ 2 d S . (4.8) 



(b) for a > 1/2, 

d (t Q J Q (t))=t a J Q _ 1 (t). (4.9) 



eft 

(c) W^e /icrae 



\J a (t)\<cr 1/2 , t>o, (4.io) 
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n 


771 


s = 1 


s = 2 


s = 3 


s = 4 


s = 5 




21 


484 


1.9324e-01 


2.6661c-01 


5.9105e-01 


1.2583e+00 


1.5732e4 


-00 


89 


1764 


8.6889e-02 


2.3262e-02 


4.7827c-02 


1.5319e-01 


2.8550e-01 


401 


6724 


4.4847e-02 


1.3919e-02 


9.1716e-02 


3.9371c-01 


6.6656e4 


-00 


1637 


26244 


1.6146e-02 


1.5521e-02 


1.4326e-01 


3.0044e+00 


7.2167e+ 


-01 



Table 6: Maximum error of MSN interpolant in region 1.1 < r. 



Proof. Part (a) is proved, except with a different notation in [HI Theorem 3.3, p. 155]. Part (b) is a 
straightforward consequence of the series expansion for J a [18, Formula (1.71.1)]: 

k=0 y ' 

The estimate (|4~TU)) follows from [HI Formula (7.31.5)]. □ 

Proof of Theorem I4.H Without loss of generality, we may assume in this proof that ||2?^/i||i.[o.oo) = 1- 
First, we prove (|4.3|) . The first equation follows immediately from the definitions and the fact that 
h(t) > for all t. Since /i(||k||/i) = if ||k|| > bt, \h{u)\ < c for u e R, and the cardinality of the set 
{k e Z' : ||k|| < bt} does not exceed c\t q , we see from the definition that < ^t(h, 0) < C2t q . This 
proves the last inequality in (|4.3I) . 

In the proof of (|4.2p . we can assume that £||x|| > 1. In this proof only, let ho(x) = /i(||x||), x € R 9 . In 
view of the Poisson summation formula |15[ p. 251] (our notation is different), we have for x G [— 7r,7r] 9 , 

x) = (2ir) q t q ho{t(x + 2nk)). (4.11) 
Let k £ Z 9 , i(x + 27rk) = y, and ||y|| = r. In view of Lemma l4~TT a). we have 

hob) = 7^72 J o H S )J {q _ 2)/2 (sr) S q ^ds. (4.12) 
Let a > —1/2. The equation (|4.9|) used with a + 1 in place of a shows that 

' J a (rs)s a+1 ds - r- Q - 2 f U J a (v)v a+1 dv = r" a " 2 (rn) a+1 J a+1 (ru) = M " +2j "+i( M ) . 



Consequently, an integration by parts in (|4.12p yields that 

/>oo 

h(s)J a {sr)s a+1 ds = r- 1 / Vh{u) J a+1 (ru)u a+2 du. 
Jo 

Repeating this Q times, we obtain 

i />oo 

h(s)J a (sr)s a+1 ds = r~ Q / V Q h(u) J a+Q (ru)u a+Q+1 du. 

Jo 

We recall that h'(u) = if u £ [a, b\. Consequently, 

\ h(s)J a {sr)s a+1 ds = r- Q / V Q h{u) J a+Q {ru)u a+Q+1 du. 

Jo J a 

In view of (|4. 10|) and the fact that a + Q + 3/2> 0, we deduce that 

/•oo 

/ h(s)J a (sr)s a+1 ds <cr- Q - 1 ^\\V Q h\\ 1 . [Q ^ ) ^cr- Q - 1 / 2 , a > -1/2. (4.13) 

.70 
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Using q/2 — 1 in place of a and substituting the resulting estimate into (|4.12[) . we obtain that 



V;(x + 2^k)) = My) < ^ 1/2 - 9/2 - Q H2^lli,[o,oo) = (f||x+ c 27rk||)jr (4-14) 

When k ^ 0, we have 

||x + 27rk|| > |x + 27rk|oo > 27r|k| OD - tt > ^k^ > 7r||k||/V«. 
Since Q > (q + l)/2, we have R> q, and hence, 

\h (t(x + 27rk))\<ct~ R E liOfl ^ c (*ll x ID" fl - (4-15) 

keZ 9 , k/0 keZ«, ||k||>l " ' 

If k = 0, yields \h {ix)\ < c(t\\x\\)- R . Together with (l4~T5|) and (|4TTT|) . this implies 

Since R > q, we see from (|4.2[) that 



/ \y t (h,x)\dx<ct q - R [ \\x\\- R dx = ct q - R [ 

J\\x\\>l/t /||xll>l/t Jll 



fOO 

u q ~ 1 u^ R du < c. 

'||x||>l/t J\\x\\>l/t Jl/t 



Since f|4.3[) shows that /i| x M <:1 / t |^t(/i, x)|dx < c as well, we have proved (14.4)) in the case when p = 1. 
The estimate (|4.4j) in the case p = oo follows from (|4.3p . The general case is obtained using the convexity 
inequality 

Hsll P <ll3ll^'llslli /p , 9e£ 1 nCK?<». 



□ 



If / G L 1 , we define 



^/.x)~7nj / /(y)* t (n,x-y)dy. (4.16) 

(27T)« 7[_ Wj7r l 3 



The following theorem summarizes some facts related to this operator. 

Theorem 4.2 Let n satisfy the conditions of Theorem \4-l\ 1 < p < oo, /eL p . 

(a) FKe Ziaue 

|k t (/l 1 /)|| p <c||I>^||i, [ o,ao}||/||p S *>Q- (4-17) 

(b) in particular, if h(t) = I on [0, 1/2] and /i(i) = on [l,oo) ; t/ien 

EtAf) < l|/-o-t(A J /)||p<c(l + ||2?^|| 1 , [ o, 00 ))^ /2 , p (/) J t>0. (4.18) 

(c) // s > and f G WJ t/ien 

£t, P (/) < ^||© Q /i||i, [0 ,oo)^,p(/ (s) ), * > (4.19) 

and 

lk t (n,/) (s) || P <c||2?«^|| li[0iCo) t s ||a t (fc,/)|| p , t>0. (4.20) 

(d) (Bernstein inequality) In particular, ift>0 and T G H', 

imiwy <rf s ||T|| p . (4.21) 

Proof. In view of gU), ||**(/i,o)||i < c||D«/i||i,[o,oo). The estimate KWi is now clear in the case 
p = oo, and follows from Fubini's theorem in the case when p = 1. An application of Riesz-Thorin 
theorem leads to the intermediate cases. 

Next, let h(t) = 1 on [0, 1/2] and h{t) = on [l,oo). Then a t (h,T) = T for all T G Hj? /2 . Therefore, 

||/ - at (/i, /)|| p = ||/ - T- <7 t (n,/ - T)|| p < c (l + ||^/i|| 1)[0)Oo) ) ||/ - T]\ p , T G H« /2 . 
This implies (|4TT5]) . 
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Next, let set, and in this proof only, g t (u) = (h(u) - h(2u))/(u 2 + l/t 2 ) 3 / 2 , u £ [0,oo), t £ (0,oo). 
Using the fact that gt(u) = if u £ [0, 1/4], it is not difficult to verify that g t satisfies the same conditions 
as ft and Hf^gtH i,[o,oo) < c\\T> ®h\\ i,[o,oo) with constant independent of t. Since at(h, f) — a 2 t(h, f) — 
t~ s ut(gt, f ), this implies that 



ht(hJ)-a 2 t(hJ)\\ p <ct- s \\f^\\ p \\V Q h\\ ltl0 ,oo), t>0, s£ 
Next, let s > 0. Then (TO!?!) leads to 



(4.22) 



EtAf) < \\f-<rt(h,f)\\p = 



'J2(a 2 k t (h, /) - <r 2 k+i t (h, /)) 



< J2 IK'tfo /) - <r 2k+H (h, f)\\ p < ct- s \\f^\\ p \\vQh\\ U0iOoh 



k=0 



If T G H< satisfies ||/ (s) - TW|| p < 2E t Af {s) ), then this implies that 

EtAf) = EtAf ~T)< ct- s WDQh\\ 1A0tOo) \\fM - T«|| p < ct- s EtAf {s) )\\V Q h\\ U0tOo) . 
This proves (gHU). 

We observe that ai(h, /) = /(0) = ai(h, / (s) ). We observe also that (|4.22[) holds also for s < 0. So, 
if s > 0, we may apply (|4.22[) with in place of / and — s in place of s to conclude that 



lk*(^/ w )-^»('»,/ w )l|p<c*'ii/y|2?^iii l[ o l0 o). 



Hence, for n > 1, 



\\a 2 J s HhJ)\\ p = \\a 2 ,(h,n 



Mil - 



CTl (ft, /«) + £ (a 2 , +1 (ft, /«) - a 2k (h, /«; 

k=0 

-1 

< \\*i(hj^)\\ p + j2\\^(hj {s) ) - ^(hj {s) : 

k=0 

I fe=o J 

This leads to (|4.20[) . The estimate (|4.21[) is obtained by using (14.201) with 2t in place of t and T in place 
of /, where we may use a fixed ft, so that the constant is independent of the function ft used in the rest 
of the statements of this theorem. □ 

Our next major goal is to prove Theorem 14.41 In this section, we develop the properties of the kernels 
fy n which are required in this proof. Let {y^f^ C [— tt, n] q , m > 1 be an integer with 

min||y, -y fc || > 1/m. (4.23) 

We note that this implies M < cm q . In the sequel, we will assume tacitly that {y^}*^ is one of the 
members of a sequence of finite subsets of [— tt, 7r] 9 . Thus, M and m are variables, and the constants are 
independent of these. If a = {afc}?*Lo ^ s an y sequence of complex numbers, we define 



a Up 



{Er=oM P } 1/P , if 1 < P < oo, 
su Pfc>oKI, ifp = oo. 



If a is in a Euclidean space R D , ||a||#> := ||(0, a\, ■ ■ ■ , a_o, 0, • 
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Proposition 4.1 Let n > 1 be an integer, 1 < p < oo, a E M. , h, Q, R be as in Theorem \4-l\ and 
G(x):=££ 
(a) We have 



G(x) :=Y!? =l a j * n (h,-x.-y j ), xG [-tt.tt]' 



||G|| P < cn^' {1 + (m/n)«} W ||2?^|| li[0!Oo) ||a|| £P . (4.24) 

(b) Suppose that there exists a compact interval I C (0, 1] and a constant cq — Co(h, I) such that h(t) > cq 
iftEl. Then there exists C\ > depending on I, cq, q, and Q such that for n > C\m, 

c 2 n-i/p'\\G\\ p < ||2?«ft|| 1)[D)00) ||a||^ < c&T*U \\G\\ V . (4.25) 

The proof requires a number of preparatory results, some of which we find of interest in their own 
right. 

Proposition 4.2 Let h satisfy the conditions of Theorem \4-l\ {yj}fLi C [— ir,Tr] q , m > 1 be an integer 
with minj^fc ||yj — y^H > 1/m. For integer n > 1 and x E [— ir,Tt] q , 

\^ n (h,^-yj)\<cn"(m/n) R \\vQh\\ h[ o, 00 y (4.26) 

J,l|x— y_j||>l/'m 

-ffence, 

M 

— ^ |\Mfc,x- yj -)| < c(n/m)« {1 + (m/n) R } ||2>«*|| li[0iOo) . (4.27) 

m 3=1 

Proof. Without loss of generality, we may assume that ||2? ( 3/i|| uo,oo) = 1- In this proof only, let 
Zfc = {j : k/m < ||x — yj\\ < (k + l)/m}, k = 1, 2, • • •. We note that since the minimal separation 
amongst y^'s does not exceed 1/m, there are at most cfc 9-1 elements in each Z/.. We note that since 
Q>(q + l)/2, R=(q- l)/2 + Q>q. In view of g^]), we have 



|* n (fc,x- yi )| < en" £ H|x-y,| 
j.llx-yil!^ 1 /™ jMl x -yj ll>!/ m 



-iJ 



= cn^^ ^ llx-yj"* < cn q - R m R J2k"- 1 - R 

k=lj£Z h k=l 

< cn q {m/n) R . 



This proves (Ej~2"o) . 

In light of (|4.23p . the number of yj's with ||x — y^ || < 1/m is bounded independently of M and m. 
Hence, (I4.3|) implies that 

^ |* n (fc,x- yi )| <cn«. 

i.llx-yill^ 1 /" 1 

Together with (Ei~2"oT) . this leads to P~2"7)l . □ 
For / : {yj} — >■ M, we will write 

i/p 



ll/l 



l p = { {^E^xl/(y,-)l p } ifi<j»< 

[ maxi< 3 < M |/(yj)|, if p = oo. 



Theorem 4.3 Lei 1 < p < oo. For any integer n > 1, and T G H^, we /lave 

|||T||| p < c(n/m)«/P {1 + (m/n) fl } 1/p ||T|| p . (4.28) 

Proof. In this proof only, let h : [0, oo) — > [0, oo) be a fixed, infinitely differentiable function, h(t) = 1 
if < t < 1/2, ft,(£) = if t > 1, and we choose Q = q + 1. The constants will depend upon this h, but 
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h being fixed in this proof, this dependence need not be specified. A comparison of Fourier coefficients 

] 



shows that for T e H|„, 



T (y) = TtAt/ / T(^ 4n (h, x - y)dy. 



In view of (I4.27[) . we obtain 

]\* in (h,X- yj )\ \ <c(n/ m y{l + (m/n) R }\\T\\ 



M ( M 

— V \T(yj)\ < \\T\\i max { — V 

TO 9^l KJJM-W » 1 x€[ _„ ]t ] m q 

3=1 t_ 3=1 



If / £ L 1 , we apply this estimate with <T2 n (/i, /) in place of T, and use Corollary 14.21 (with p = 1) to 
deduce that 

Hkan(/i,/)|||i < c(n/m)« {l + {m/n) R } \\f\\ x . 
In view of (|4.17[) (with p = oo), it is clear that for / G 

Hka»(/»,/)|||oo < ||<7-a»(/i,/)||oo < cll/Iloo. 
An application of Riesz-Thorin interpolation theorem now implies that for 1 < p < oo, and / G L p , 

\\W2n(h,f)\\\ P < c(n/m)^ {1 + (m/n) R } 1/p \\f\\ p . (4.29) 

If T G H«, then cr 2 „(/i, T) = T. Therefore, (TO?)) implies fl3~J5J}. □ 
Proposition 14. 31 below is perhaps well known. A proof can be found in [71 Proposition 6.1]. 

Proposition 4.3 Let M > 1 be an integer, A be an M x M matrix whose (i,j)-th entry is Aij. 
1 < P < oo, and a € [0, 1). If 

M M 

< a\A jd \, Y, I Ail < a\A jd \, j = l,---,M, (4.30) 

i=l i=l 

and A = mini<i<M > 0, then A is invertible, and 

l|A _1 b|| £ p < ((1 -a)\y l \\h\\ iP , beR M . (4.31) 



We are now in a position to prove Proposition 14.11 
Proof of Proposition 14.11 Without loss of generality, we may assume that ||2? Q /i||i.[o.oo) — 1- In 
view of (|4.27[) . we have for x e [— ir,Ti] q , 

M M 

|G(x)|<2|a J ||* n (A,x-y i )| < ||a||<~ |*„(fc,x - yj )| < cn« {l + (m/n) R } ||a||^. 

3=1 3=1 

Thus, 

||G|U <cn 9 {l + (m/n)*}||a||*ao. 

Using (|4.4[) . we see that 

M 

||G||i <X)KII|*n(^°-yj)||i < c||a|| £ i. 

3 = 1 

An application of Riesz-Thorin interpolation theorem with the operator a 1— > Y]j—i Q"i^n(h, ° — y/) 
implies (|4~24|) . 

Next, if the hypothesis in part (b) is satined, then 

*»(M)> Y H\m/n)>cni. (4.32) 

k, ||k||/neJ 
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Therefore, flUH} and (|4~26| show that for n > Cirri, I = 1, • • • , M, 

Y,\^n(h,y e -y 3 )\ < (1/2)*„(M). (4.33) 



3 = 1 



In this proof only, let A be the matrix whose (£, j)-th entry is ^f n (h,yi — yj) and b e K M be defined 
by be = G(y^), ^ = 1, • ■ ■ , M. In view of (|4.33l) . (14.301) is satisfied with 1/2 in place of a, and in view of 
(|4.32[) . we may choose A to be cn q . Hence, Proposition 14.31 implies that A is invertible, and 

HA-^H/p < cn- 9 ||b||^. 

Since, A _1 b = a, we have proved that 

||a||ip < cn- 9 ||b|| fr . (4.34) 
Since G £ H*, we obtain from Theorem 14.31 that 

|||G||| P = m-*S*\\b\\t, < c(n/m)«/" {l + (m/n) R } 1/p \\G\\ p . 
Since n > C\m, this gives 

\\h\\» < cn q/p \\G\\ P . 

Together with (|4. 34[> . this leads to the second inequality in (I4.25p . The first inequality follows from (|4.24l) 
and the fact that n > C\m. □ 



4.2 Sobolev kernel 

Our goal in this section is to prove Proposition 12.11 and Theorem 14.41 and establish a few other facts 
regarding the kernel K s . In particular, we will give in Theorem 14.51 an estimate for the norm of the 
interpolation matrix K2 S (yj — yk)- In the sequel, we assume Q > (q+l)/2 is an integer, h : [0, oo) — > [0, oo) 
is a fixed, Q — 1 times continuously differentiable function with an absolutely continuous derivative 
frW -1 ), h(t) = 1 if < t < 1/2, h(t) = if t > 1, and h is nondecreasing on [0,oo). We will write 
g{t) = h(t) — h(2t). Since h is fixed, the dependence of various constants on h need not be indicated. For 
s£l, we will write 

* n ,.(x) := J2 2(l|k||/2")(||k|| 2 + I)'*' 2 exp(ik • x). (4.35) 

kez<! 

The following lemma lists some interesting properties of ^ n ,s- 
Lemma 4.2 Let s£i. We have 

2n(q-s) 

|Vl/ - (x)l ^ min(l, ( 2»||x||n ' * ^ 



Further, 



and for 1 < p < oo , 



max |#„, s (x)| = * n , s (0)~2"(«- s ), (4.37) 

X£ [— 7T,7r] ? 

||*„,,|| P <c2"^'- s ). (4.38) 

PROOF. In this proof only, let g n (t) = g{t)/(t 2 + l/n 2 ) s / 2 . Then for x e [-tt, 7r] 9 , 

*„,,(x) =2-" s * 2 „(.g 2 , 1 ,x). (4.39) 

Each g n satisfies the conditions of Theorem 14.11 with a = 1/4, 6=1. Moreover, ||2?^<7nlli,[o,«0 — c - 

Therefore, all assertions of the lemma, except for the second relation in (|4.37[) , follow directly from 
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TheoremEU Theorem O also implies that #„, s (0) < c2 n ^~ s \ Since ,g(l/2) = h(l/2) - ft(l) = 1, and 
g is continuous, there exists a nondegenerate interval / C [1/4, 1] such that g(t) > 1/2 if t € I. Hence, 

*„,.(0) - 5(l|k||/2")(l + ||k|| 2 )- s / 2 > 5(l|k||/2")(l + ||k|| 2 )- s/2 > c2 n ^~ s \ 

kezi kezi, ||k||/2™e/ 

This proves the second relation in (|4.37l) . □ 
Proof of Proposition 12.11 Since s > q/p, (|4.38[) used with p' in place of p shows that 

oo oo 

<cJ22 n(q/p - s) <oo. 

n=0 n=0 

So, the sequence of trigonometric polynomials, defined by 

N N 

P N (x) = 1 + J2 = 1 + E E 5(l|k||/2")(l + ||k|| 2 )- s / 2 exp(*k • x) 

n=0 n=0keZi 

converges in L p . All the sums in the above expression being finite sums, we obtain for N > 0, 

N 

P N (x) = l+Y,J2 3(Hk||/2 n )(l + ||k|| 2 )^/ 2 exp(zk ■ x) = £ h(\\k\\/2 N )(l + ||k|| 2 )-/ 2 cxp(zk • x). 
kez-j n=o kez q 

If k e Z q , and 2 N > 2||k||, then ftfllkH^) = 1, and P N (k) = (1 + ||k|| 2 )- s / 2 . Denoting the L p ' -limiting- 
function of Pjv by K s , it follows that K s <E LP and satisfies (12.41) . Moreover, Pjv = c 2 iv (ft, if s ), and the 
bound on H^n.sllp' m fl4.38|) used with p' in place of p shows that 

\\K s -a 2 N(h,K s )\\ pl <c2 N ^- s \ # = 0,1,2,..-. (4.40) 

Both sides of the first equation in (|2.5[) have the same Fourier coefficients, and hence, they are equal 
almost everywhere. Similarly, a comparison of Fourier coefficients shows that if s (— x) = if s (x) for almost 
all x. This implies the second equation in (|2.5|) . 

For / £ Wf , a comparison of Fourier coefficients again shows that for integer m > 0, 

a 2 m (ft, /, x) = — L- / a 2m (ft, If,, x - y)/( s ) (y)dy. 

So, (|QU|) implies that 

Hence, the series a\{h, /) + J2m=i( tT ^ m f) ~ cr 2 ™-i (h, /)) converges uniformly. It is clear that this limit 
is almost everywhere equal to /, and by choosing the continuous representer in the equivalence class of 
/ to be /, the limit is /. Moreover, 

oo 

Sa»,oo(/) < 11/ (Kf) Hoc < E h^(hJ)-^^(h,f)\\oo < 2™(^- s )||/W|| p . 



This implies the first estimate in (|2.6I) is now clear. The second set of estimates are proved similarly. □ 

Our proof of Theorem l2. II requires the following theorem that describes an approximation of a typical 
element of the span of {K s {o — y,)}. We recall that the solution of the minimization problem (12.7[) is in 
this span (with 2s in place of s). 

Theorem 4.4 Let 1 < p < oo, s > q/p, {aj}f =1 C M, G(x) = Y.f=i ajK s (x - Yj), x € [-7r,7r]«, and 
m > 1 be the smallest integer such that miny^.^y^ ||yj — yfe|| > 1/m. Then there exists an integer N* , 
independent of G, such that N* ~ m and 

\\G-a N *(h,G)\\ p , < (l/2)||G|| p ,. (4.41) 
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Proof. As in the proof of Lemma[OJ in this proof only, we write g n {t) = g(t)/(t + l/n) s . Then each 
g n satisfies the conditions of Theorem 14. 11 with a = 1/4, 6=1. Moreover, \\T> ^ g n \\i,[o,oo) ~ 1; an d (|4.39l) 
holds. In this proof only, let 

M 

G„(x) := y^a 3 § W|S (x - yj) = a 2 « (h, G, x) - <7 2 n-i(/i, G,x), x e [-tt, tt] 9 . 
Then (|4.17l) implies that ||G„|| P < < c||G|| p '. Moreover, the proof of Proposition I2TT1 shows that 

oo 

G(x) - a 2 N {h, G, x) = G «( x )' ( 4 - 42 ) 

n=Af 

with convergence in the sense of L p . 

In view of (|4.39[) . (|4.25[) applied with ^ 2 ^(g n ) yields that for n > log 2 (Gim), 

C2 2^-i/p)\\G n \\ p , < \\a\\ iP , < c2"( s -^}||G„|| p ,. (4.43) 

We now choose L so that 2 L is the smallest power of 2 exceeding C\m. Then the second inequality in 
(|4.43p . used with L in place of n, gives 

||a|| £p < < cm<< s ~^\\G L \\ p , < cm^-^WGWp,. (4.44) 

From (|Q2]1 . fl~4Til) . and (|4~44l) . we conclude that for 2 N > Gim, 

OO OO 

||G-a 2 *(/i,G)|Ly < \\Gn\\ P ' < c\\ a \\ f y £ 2-"^^ < c^-^^-^HGH^. 

We now choose TV so that 2 N ^ m and the last term above is at most (1/2)||G|| P >, and set TV* = 2 . □ 

We note a consequence of the proof, which might be of independent interest in view of the fact 
that the interpolant which yields the minimal Sobolev norm amongst all interpolants is in the span of 
{K 2s (o — Yk)}kLi- The following theorem gives the norm of the inverse of the interpolation matrix 
{K 2s {yj — Vk)) in terms of the minimal separation (equivalently, m). 

Theorem 4.5 Let s > q/2, and I be the M x M matrix whose (j, k)-th entry is K2 S {yj — yk), where the 
points yj satisfy fa4-£>ty - Then X is positive definite, and 

< cm 3 - 1 ' 2 . (4.45) 

Proof. We observe that a comparison of Fourier coefficients shows that 

K 2s {y 3 - yk) = 77t~tt / K s (x - yj)K s (x - y k )dy. 
(2tt)« J [ _ 7ri7r]g 

Let a e R M , and G = Ylj=i ajK s (o — y rj ). Then the above identity leads to 



3- 

2 



M 1 

a J a k K 2s {y J - y fe ) = —- ■ \ 
k,j=l V ' ' 



M 

J2 a 3 K s(*-yj) 



dy = \\G\\i (4.46) 



The estimate (|4.44|) used with p = 2 now shows that 

M 

J2 ^a k K 2s { yj -y k ) > cm-^/ 2 >||a|| £2 . 

k,j=l 

Thus, I is a positive definite matrix. In view of the Raleigl^Ritz theorem [?J Theorem 4.2.2, p. 176], the 
lowest eigenvalue of this matrix is at least cm - '^^ 2 '. This implies (|4.45l) . □ 

Although not strictly a property of the kernels K s , we find it convenient to record the following lemma, 
which will be needed in our proof of Theorem 12.21 This lemma is proved in much greater generality in 
[5j Theorem 3.2, Chapter 15]. 

Lemma 4.3 Let 1 < p < oo, s' > 0. Then for any c > 0, the sei 73 C;S <, P := {/ 6 L v : sup n>1 2 ns E n ^ p (f) < 
c} is compact in L p . 
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4.3 Background on approximation theory 

The proof of Theorem 12.21 depends upon a number of facts from classical approximation theory, as well 
as some recent developments. In this section, we review the necessary facts. 

First, for integer r > 1, the modulus of smoothness u> r (f,5) of a 27r-periodic univariate function 
/ £ L p {{— 7r, 7t]) is defined by first defining the forward difference operator 

r 



k=0 

and setting 

w r (f,S) := max||AJ7|| p . 

\t\<s 

If / £ L v {[— 7T, 7r] 9 ) and r = (ri, ■ ■ • , r g ) > 0, r ^ (0, •••,0), is a multi-integer, then the modulus of 
smoothness is defined in [T9l Section 3.4.34] by 



u) T (f. h) = max 

|ti|<fci,-,|t«l<fc« 



n *z* f 



\k=l 



where the notation A£* k means that the operator A£* is applied to the fc-th variable in the argument 
of / and A° fc k f means that no difference is taken with respect to th fc-th variable. We will write ejt to 
denote the vector in R 9 with k-th coordinate equal to 1 and the remaining coordinates equal to 0. 

For an integer n > 0, the class of (rectangular) trigonometric polynomials of order at most n is defined 

by 

W q n R := | a k exp(ik- (o)) : a k € 

For f E L p , the degree of approximation from Mf,' R is defined by 

E n , P ; R (f)= inf J|/-P|| P . 

Let r > be an integer, a £ (0, 1] .It is proved in [TJ2 Section 3.6.4] that for / £ L p , the relation 

supn r+a E ntP . !R (f)<oo (4.47) 

n>l 

holds if and only if / has almost everywhere defined partial derivatives D r k f satisfying 

max supS^ a uj2 ek (D r k f,Se k ) < oo. (4.48) 

l<fe<9 «5>0 

(The formulation in [19] is not quite precise. However, the version which we have stated can be obtained 
using the same ideas as in [T5]. See [501 4.5.6] for an analogous statement in a slightly different context.) 
Since H« C C H^ n , we have 

E^qn,p{f) < En,p;R(f) < E n ^ p (f). 

Thus, (|4.47l) is equivalent to 

supn r+Q i?„.p(/) < oo. 

n>l 

We summarize these observations in the following proposition. 

Proposition 4.4 Le£ 1 < p < oo, and f £ L p , r > be integer, and a £ (0, 1]. Then 

sup n r+a E n , p (f) < oo (4.49) 
n>l 



if and only if f has almost everywhere defined partial derivatives D k f satisfying {J^J+Sty . 
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Next, we recall some results from the theory of algebraic polynomial approximation. Let II 9 , denote the 
set of all algebraic polynomials of coordinatewise degree at most r. We wish to construct an approximation 
to a continuous function / on [— 1, l] 9 , defined analogously to (|2.1I) . based on an arbitrary data set 
C C [-1, l] 9 . The mesh norm S c := 5(C, [-1, l] 9 ) is defined analogously to (JH]). We divide [-1, 1]« into 
equal subcubes of side 25c i the set of these subcubes will be denoted, in this part of the discussion only, 
by IZc- Each of the subcubes has at least one point of C. We form a subset C\ C C by choosing exactly 
one point £ in each £E TZc- Then it is clear that Sc t ~ Sc- In the following discussion, the points in 
C \ Ci do not play any role, and accordingly, we rename C\ to be C. The following proposition follows 
from [9j Theorem 3.1], by taking the functional P — > L ^ q P(x)<ix in place of 7 in that theorem. 

Proposition 4.5 Let C, TZc be as above, r > 1 be an integer. There exists a constant 7 := 7(7-, q) with 
the following property. If 5q < 7, then 

£vol 9 (i^)|P(e)|~ / |P(x)|dx, P6H1. (4.50) 
Zee J [~W 

Further, there exist real numbers {a^ : £ G C}, such that 

K| < cvol 9 (P 5 ) < cS*, £eC, (4.51) 

and 

E a ? F ^)=/ F « dx > ^en^. (4.52) 

tec J l-^ q 

In this section only, let pk denote the orthonormalized Chebyshev polynomial of degree fc, with positive 
leading coefficient. We define 

v r(x,y) = ^2pk(x)pk(y), x,y e [-1, 1], r = 1,2,- • • . 
and extend this definition by writing 

v r (x,y) = Y[v r (x e ,y e ), x,ye[-l,l] 9 , 

i=\ 

If / : [—1, l] 9 — 7- R is continuous, C and {a^} are as in Proposition 14. 5[ we define 

K-(/,x) = ^o c /(0wr(x,0- (4-53) 

Using (|4.51l) . (|4.50j) . we conclude that 

||W)||oo,[-i,i]« < cWII/IU.i-n].. (4.54) 

In view of g3U), U r (P) = P for all P e IT 9 . So, choosing P* G II 9 with ||/-P*|U,[-i,i] a < 2inf Pen? ||/- 
P\\oo, 1431 yields 

11/ - W)||oo,[-l,l]« = 11/ - P * ~ Mf - )||oo,[-l,l]« < C||/- P*||oo,[-l,l]« < C ^f ? ||/ - P||oo,[-l,l]«- 

(4.55) 

Using the direct theorem of approximation theory [191 Section 5.3.1], we conclude that if / has continuous 
partial derivatives of order up to r, then 

11/ - W)||oo,[-i,i]« <c(r) J2 W2e»(0£Wr), (4.56) 

fe=i 
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where the modulus of smoothness is defined analogously to (|4.48[) , except that the maximum is taken for 
only those values oft\,---,tk which don't take the argument out of the cube in question. We note again 
that the operator V r is determined entirely by the values {/(£)}{eC- 

We end this section by recording another observation, establishing a connection between the discrete 
norm used in the statement of the minimization problem (|2.1ip and the continuous L p norm, which will 
be needed in the proof of Theorem 12. lf b). 

Lemma 4.4 For integer n > 1, 1 < p < oo, and T g we have 

E |T(27rk/(3n))|4 ~\\T\\ P . (4.57) 

[ 0<k<3n-l J 

Proof. When q = 1, (|4.57|) is the classical Marcinkiewicz-Zygmund inequality [2TJ Chapter X, Theo- 
rems 7.5, 7.28]. If T G then T e H(„.... !Tl ). So, in the case when q > 1, one obtains (|4.57|) by applying 
its univariate version to each of the variables separately. □ 



5 Proofs of the main results in Section [21 

Our proof of Theorem 12. lf a) relies upon the next lemma, proved in |8| Theorem 2.1]. 

Lemma 5.1 Let X be a normed linear space, V C X be a finite dimensional subspace of X , X* be 
the dual space of X, {Xj}jL 1 C X* , and Z* be the span of {x*}^^. Suppose that the restriction map 
S : z* 6 Z* i— > z* \v is injective, and the operator norm 1 1 ^S' 1 1 1 < k for some k > 0. Then for every f £ X 
and K\ > k, there exists T{f) G V such that 

z*(T(/)) = z*(f) for every z* E Z„ (5.1) 

and 

||/-T(/)||x<(l + Ki)mf \\f-v\\x. (5.2) 

v£ V 

We will use this lemma with Wf in place of X, M 9 N , in place of V for a suitable N* , and point 
evaluation functionals in place of x* 's. 

PROOF of Theorem 12.11 The proof of this theorem is similar to that of Theorem 3.1], except that 
the details are much more complicated, requiring the use of Theorem 14. 21 and Theorem 14.41 In this proof 
only, let X = Wf, x*{f) = f(yj), j = 1, • • • , M. Since s > q/p, Proposition 12.11 implies that x* e X* , 
j = 1, ■ • • , M. Let {a 3 }f =1 C K and z* = Y%=\ a 3 x *y Let N * be as in Theorem EH and V = M q N ,. To 
estimate US'" 1 !! for the operator S as in Lemma f5.1[ we need to find T € H^, for a suitable N* , and 
estimate \z* (T)\/\\T\\ w p from below. Let /* be chosen so that ||z*||x- < (4/3)|«*(/*)| and \\f*\\ w p = 1. 
We will prove that a N ,\h,f*) e (cf. (|4~T6| ) satisfies 

\z*(T)\ \z*(a N «(hJ*))\ 
p TFrTi — - iu (h f*\\\ — ^ C II Z H**- ( 5 - 3 ) 

T£V IM || wf \\VN'{n-,J )\\wi 



The part (a) of the theorem will then follow from Lemma 15.1 
We start by observing that for / e Wf, (|2.5|) shows that 

z ' {f> "WysL,, { 5>*'fr-">} I G w/ <s, w*. 



■7T,7r ]o 



where G is defined as in Theorem 14.41 In light of the duality principle and the definition (|2 .3[) . we see 
that 

=sup{|z*(/)| : \\f\\ w? =l} = \\G\\ p ,. (5.4) 
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Further a comparison of Fourier coefficients implies that for any integer n, 

z*(a n (h,f)) = -L- / G (y)(* n (hJ)f\y)dy = -±- f G(y)a n (h, /«, y)dy 



1 



a„(/i,G,y)/( s )(y)dy 



Then 



(G(y)-a Jv .(/ 1 ,G,y))r (s) (y)dy 



(27T)9 

< HG-^.^G)^ < (l/2)||G|| p , = (1/2)11*1*. < (2/3)|«*(D|. 

Therefore, 

|*>Wfc.r))l > (l/3)k*(DI > (V4)lk*|U.. (5.5) 
Moreover, (|4.17p implies that 

hN*(h,n\\wf = \\Mh,nn P = ik^^rwjiip < c|ir«n P = c . (5.6) 



The estimate (jOjl follows from ([53]) and (|5T6|) . 

We note that necessarily, ||P(/)||wj < ||/||w?- Therefore, part (b) is a simple consequence of 
Lemma 14.41 □ 

Proof of Theorem 12.21 

To prove part (a), we observe that in view of (12.6[) and the fact that ||P*|| w p < c||/||^p for all n, the 
sequence {PJj} C B c s _ q / p oo for a suitable constant c. Let Ai be any subsequence of A. Then Lemma l4.3l 
shows that the sequence {P*} ne Ai has a subsequence {P£}neA 2 j which converges uniformly. Let P be 
the limit of this subsequence. We will show that if (|2.12j) is satisfied, then P(xo) = /(xo). Let e > be 
arbitrary. Since P and / are continuous on [— 7r, 7r] 9 , there is 5 > such that 

|/(x)-/(y)|<e/3, |P(x)-P(y)|<e/3, for all x, y G [-tt, tt]«, ||x-y||<5. 

Further, there exists N so that n > N, n G A 2 imply that ||P - P* ||oo < e/3. In view of ([2TT2]) . there 
exists n G A 2 , n > N such that some point yj jW G Y" satisfies Hy^.n — xq|| < S. Then f(yj >n ) = PJUyj,".); 
and we have 

|/(x ) - P(xo)| < |/(x) - f(y j<n )\ + \}{y hn ) - K{y 3 , n )\ + \Kbi») ~ P<Xi,n)\ + l p (y.>) - ^( x o)| 
< e/3 + + ||P,; -P||oo+e/3< e. 

Since this is true for every subsequential limit of P* , n G A, this proves part (a). 

To prove part (b), let r be an integer and a G (0, 1] be chosen so that s — q/p = r + a. Since 
PJj,/ G B CiS _ q / Pt00> Proposition 14.41 implies that they both have r derivatives satisfying (|4.48|l . In this 
proof only, let P(y) = P* (x + 6y), f(y) = /(x + Sy), y G [— 1, l] q . Then the assumptions of part (b) 
ensure that we can construct the operator V r as in (|4.53l) based on Y n n K in place of C. In this proof 
only, if k G {1, • • • , q}, F = D r k f, then D r J{y) = <5 r P(x + Sy). Further, 

A? /r>fc (P>£/)(y) = (A? /rifc F)(xo + <5y). 

Hence, 

ll^kPfc/JIUt-M]. = IIA^fcFCxo+o)!!^^, <u) 2 (Dy,(5/r)e k ) < cS a . 

A similar estimate holdes also for P in place of /. Using (|4.56p and the fact that r + a = s — q/p, we 
deduce that 

II/- K-(/)IL,[-i,i]« < cS s '^, ||P- V r (P)||oo,x < c5 s "^. 

We now observe that V r (f) = V r (V), and hence, the above inequalities imply that ||/ — P||oo,[-i,i]9 < 
cS s - q / p . Scaling back to the original scale, we obtain (|2.13[) . □ 
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